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Abstract 

A phantom network of Gaussian chains far from the point of gelation can 
be described as a gas of interacting particles represented by the cross links. 
The type of particles varies with the network functionality, whereas the type 
of interaction depends on the properties of the connecting chains. In a mean 
field approximation the Equation of State can be calculated using Mayer's 
cluster expansion. The resulting isothermal compressibility is compared for 
different cross link functionalities. 
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I. INTRODUCTION 

Even nowadays the statistical mechanics of a polymer network is far from being fully 
understood. Usually a Gaussian network is described by a model proposed by Edwards in 
1968 [|l|]. In his approach the cross link constraint leads to an effective short range interaction 
term between free polymer chains. This purpose is achieved by summing over all possible 
cross link configurations before performing the sum of configuration of the junctions. In 
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fact this model is only reasonable if used in the context of quenched cross links. Using the 
replica trick to determine the system's free energy, the calculation turns out to be tedious 
since the cross link-induced action term couples all replicated configuration spaces [Q] . 

Still there is an approach used by H. Brodowsky and S. Prager in 1963 [[J, who calculated 
the statistical mechanics of gel formation in an ideal phantom network. The phantom 
network is a model introduced by H. James and E. Guth in 1943 fU, which describes the 
phenomenological behaviour of a network, while assuming that the junctions may fluctuate 
freely without being restricted by topological constraints due to the existence of neighbouring 
chains. Prager and Brodowsky are using a microscopic ansatz, where the network is dealt 
with as a gas of free cross links interacting by means of the connecting polymer chains. It is a 
phantom network, since the chains have no excluded volume term. Thus different monomers 
can be placed on the same site. The chain configuration average has been performed before 
averaging the cross link configuration and the influence of the junctions is absorbed in an 
effective potential between cross links. In this context the cross links are treated as being 
annealed. I.e. the time of spliting up and rearranging a cross link is small compared to 
typical time scales of observation. Due to the annealed limit, topological constraints caused 
by entanglement of network chains are disregarded. 

II. CALCULATION OF THE EQUATION OF STATE 

In this article we demonstrate how the Brodowsky and Prager model changes if the 
annealed phantom network is composed of cross links having a predefined functionality. 
For clarity we assume that all connecting chains in the network are of equal length. Re- 
turning to the general case of varying chain length presented by Brodowsky and Prager is 
straightforward. 

The network consists of M cross links and N bifunctional monomer units forming K 
active chains or junctions that are attached on both ends to different cross links. Since the 
length of the junctions is fixed at n monomer units the total number of monomers is given by 
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N = nK. If only active chains are taken into account the configurational partition function 
reads 
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The function q is the number of configuration that a junction a connecting cross link i to 
cross link j can assume . For a Gaussian chain it is given by 
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This is the number of different configuration of a chain of n units whose ends are fixed at 
points a distance pij =\ — | apart. The elastic properties of this chain are determined 
by the Kuhn length Zq. 

The summation extends over all ways of attaching the chains to the M available 
cross links, a is a matrix with giving the number of junctions connecting the cross link % 
to the cross link j. Conservation of the total number of active chains is taken into account 
by the first constraint 5y^ K . 

The last term defines the number of chains starting from a cross link i. It is given by 
the functionality /. This is assumed to be an even number to assure the connectivity of the 
network. To express the constraint in a simplified way two possible cases which can arise are 
distinguished. If / is taken to be a single number, the Kronecker delta of the last constraint 
in Eq.(|I|) can be written using Grassmann variables rj that are defined by 

dt] r drjrj 
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These properties allow us to rewrite the Kronecker delta as 
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If the functionality is abitrary but even / = 2k [], the Kronecker delta can easily be expressed 
using classical spins s» = ±| 

AT M 
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where the sum extends over all possible spin configurations. For clarity the latter case 

W 

will be used for our calculation that is performed following |J. In a first step the partition 
function 
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is determined. The chemical potential v is associated to the number of active chains K. 
The summation over can be performed and yields an exponential function. Thus the 
partition function reads like the partition function of a gas 

M 
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composed of charged particles. These are interacting via an inverted Gaussian potential 



where the typical length scale is given by the gyration radius of the junctions R g = ^3nll/2. 
In the more detailed case of a fixed functionality where the classical spins are replaced by 
Grassmann variables, the network is mapped on a gas of particles that are composed of / 
fermions. In the case of variable length of junctions the potential changes to a Debye-Hiickel 
kind of interaction ||. 

The partition function of Eq. ([?]) cannot be calculated rigorously but a virial expansion of 
the Equation of State can be obtained using Mayer's cluster expansion |J. It is an expansion 
of the system's grand partition function 



X NB. Like in the work done by Brodowsky and Prager the cross links are never saturated. But 
contrary to them we assume that the functionality has to be even at any cross link. 
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in powers of gij = (1 — exp(— (3siSj<&(pij)). Thus the Equation of State can be written in a 



series in powers of the average volume per cross link v — 44 
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To determine the Equation of State up to the second virial coefficient the following 
irreducible cluster integral has to be determined 
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where r = r/R g . The value of the first cluster integrals for the case of abitrary but even 
functionality is given by 
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The integral can only be evaluated by means of numerical methods. Using the description 
for a fixed functionality /, the first cluster integral simplifies due to the reduction of the 
exponential function to a single Taylor coefficient and reads 
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The summation over all functionalities, starting with / = 2, leads to the general result in 
Eq. (|I2"D . It has to be pointed out that the relation /3\ < holds in both expressions of the 
cluster integral flT2"D and (JT3|). 



III. THE VOLUME COMPRESSION MODULUS 



Having determined the Equation of State, the volume compressibility of the phantom 
network can be calculated. One has to keep in mind that the annealed phantom network 
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is a thermalized gas. Therefore the elastic constants that are connected to an anisotropic 
deformation such as a shear modulus do not exist. Using the Equation of State flTUD the 
isothermal compressibility i.e. the relative change of volume under isotropic compression 

can be calculated. Expanded up to the second virial coefficient the inverse compressibility 
reads 

^ 4-^(5)" (15) 
That allows to compare three different regimes of the network. 

I. The gyration radius may be much larger then the typical distance between cross links, 
given by v ~ 3 . That is a regime where cross links far apart are connected by a junction, 
Rg 3> v. Thus the connecting chains are highly entangled so that instead of growing 
linear with the volume, is now growing with a quadratic dependence on the volume 
V 2 . Still, due to the highly entangled junctions, there should be a correction by the 
intrinsic topology of the network, that we expect to be of the same or of higher order 
than the phantom network's contribution. 

II. If the gyration radius is much smaller than the cross link distance Rg <C v, the network 
mainly consists of small subnetworks, constituted of two cross links and / junctions. 
Since the density of these subnetworks is very small - their extension is also given by 
R g - the system behaves like an ideal gas and fin? 1 — -■ 

III. It is essential to point out, that the regime, where the assumption of a phantom 
network should give the main features of the networks compressibility is within the 
region R? g ta v . The inverse compressibility is now given by an ideal g with a 

rescaled temperature, that strongly depends on the functionality of the network (cf. 
Eq. 
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These results are not comparable to the elastic equation of state, obtained for phantom 
networks using a Flory approach 0. Since the anisotropic deformation cannot be treated 
by mapping the network on a gas, only homogeneous volume change are described by the 
model presented in this paper. The dependence on the functionality of the networks in the 
first virial correction differs from the results of the linear compressibility modulus due to 
anisotropic stress of the phantom network as given in j7|. 

Experimental or numerical results to find the volume compressibility of a phantom net- 
work are not known to us, so that we could only speculate about its contribution to a real 
polymer network. In a real system, there should also be important contributions given by 
the soft- (or hard-) core potentials of the molecules, as well as contributions due to en- 
tanglement effects. Still, the model of Brodowsky and Prager allows to find the volume 
compressibility using a very general and straightforward scheme. 

IV. ACKNOWLEDGMENTS 

This work has been partially supported by the Deutsche Forschungsgemeinschaft schu 
934/1-3, NSF Grant DMR 962 8224 and the Donors of the Petroleum Fund of the American 
Chemical Society. 



REFERENCES 
[1] S.F. Edwards, J. Phys. A 1, 15 (1968). 

[2] P. Goldbart and N. Goldenfeld, Phys. Rev. A 39, 1412 (1989). 

[3] H. Brodowsky and S. Prager, J. Chem. Phys. 39, 1103 (1963). 

[4] H.M. James and E. Guth, J. Chem. Phys.ll, 1733, (1943). 

[5] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, Oxford, 1993. 

[6] J. Mayer, J. Chem. Phys. 5, 67 (1937). 

[7] J.E. Mark and B. Erman, Rubberlike Elasticity, A Molecular Primer, New York, 1988. 



8 



